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We present a general formulation to suppress pure dephasing by multipulse control. The formula 
is free from a specific form of interaction and is expressed in terms of the correlation function 
of arbitrary system-reservoir interaction. We first apply the formula to a phenomenological two- 
level model where the correlation function of the interaction decays exponentially. In this case, we 
analytically show that the pure dephasing time is effectively lengthened by the multipulse control. 
Secondly, we apply the formula to the spin-boson model where a spin nonlinearly interacts with a 
boson reservoir. We find the multipulse control works well when the pulse interval is sufficiently 
shorter than the correlation time of system-reservoir interaction. Moreover, in this case, the pure 
dephasing can also be suppressed by adjusting the pulse interval to the period of dynamical motion 
of reservoir. 

PACS numbers: 03.65.Yz,03.67.Hk,05.30-d 



I. INTRODUCTION 

The reliability of quantum information rests on the sta- 
bility of quantum superposed state. Since the quantum 
superposition is vulnerable to the effect of the environ- 
ment, the degradation of the superposed state would af- 
fect the accuracy of the quantum information, especially 
in quantum computation. 

In order to reduce errors caused by the environmental 
effects, many methods have been proposed. These are 
roughly divided into two categories whether the meth- 
ods requires ancillary bits or not. Since the first pro- 
posal by Shorjl| , considerable works on the methods with 
additional qubits have been done to overcome the error 
by distributing the information of a qubit over several 
qubits [| , js], flj , ||, ||, including experimental studies of 
NMR[g~9|, flO| and trapped ions jd). However, in these 
methods, it has been a main problem to keep the coher- 
ence of all the qubits including ancillary bits, which may 
be an obstacle for scale up of quantum computer. 

Two distinctive approaches have been taken to attain 
strategies for control of decoherence without additional 
bits. Firstly, closed-loop control (quantum feedback) 
method has been proposed |l2|, [ll| [l4|, [l5| . Since they 
require the measurement on the system to control the 
error, the quantum detection efficiency gives a limitation 
to execution of the control. The other method to sup- 
press the error by decoherence without ancillary bit has 



[P?],[Tf|. In the pio- 



been called the open-loop control |1 6, 
neering works of Viola and Lloyd and Ban|17[, they 
have proposed that a pulse train can be used to elimi- 
nate pure dephasing of a spin caused by interaction with 
a quantum reservoir. They assume that the spin linearly 
interacts with boson reservoir. The pure dephasing of 
superposed state is successfully suppressed by applying 
many successive short 7r pulses. An actual evaluation has 



been done by using the coupling function for the macro- 
scopic tunneling phenomena as proposed by Caldeira and 
Leggett[^9). Although the pulse methods have been well 
known in the field of NMR to decouple unnecessary spin- 
spin interaction |^|, their formulation can also be applied 
to suppress decoherence in various systems other than the 
nuclear spin system. The applicability has been gener- 
ally discussed in Q . Application of a 2ir pulse train to 
a two-level system, which interacts with a boson reser- 
voir by continuous coupling function, has been studied 
to discuss the controllability of unwanted spontaneous 
emission [^2| |23fl . The magnetic state decoherence by col- 
lisions in a vapor has been analyzed to be suppressed by 
ultra fast pulses^, ^5|. Moreover, the method is applied 
to the damped vibrational mode of a chain of trapped 
ionsQ. 

It is however noted that, in the analysis |l6|, [l7|, the 
interaction between the spin and bosons are assumed to 
be linear in the boson field amplitude. We should re- 
call that in many systems the dephasing phenomena can- 
not be described with the linearly interacting spin-boson 
model, except for the ohmic case of the Caldeira-Leggett 
model. This is because the dephasing rate is given by 
the zero-frequency component of the power spectrum of 
interaction amplitude |27], 5s|. A typical example is a 
localized-electron phonon system in which the phonon 
with zero frequency (i.e., translation of a whole crystal) 
does not interact with the ^ spin (the two-level electron) . 
The nonlinear interaction between the spin and phonons, 
where simultaneous multiple absorption and emission of 
phonons occurs, should be taken into account for the cor- 
rect description of the dephasing phenomenon. In order 
to discuss the effectiveness and generality of the multi- 
pulse control of qubits, a spin-boson model that includes 
nonlinear interaction is necessary. 

In this paper, we present a general formulation free 
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from a specific form of interaction, and apply this to the 
two cases (l)a phenomcnological two- level model where 
the correlation function of the interaction decays ex- 
ponentially and (2)a quadratically interacting localized- 
electron phonon system. 

This paper is composed as follows: In section II , a gen- 
eral formula which describes spin dynamics under mul- 
tipulse control is obtained for an arbitrary interaction 
with the cumulant expansion method. The expression is 
written with the correlation function of interaction am- 
plitude. In section HI, we first apply the formula to 



a simple case of the exponentially decaying correlation 
function in order to generally discuss the memory effect 
of the reservoir causing the spin dcphasing. We next ap- 
ply the formula to a quadratically interacting localized- 
electron phonon system, and the dynamical behavior of 
spin is analyzed. Discussion and concluding remarks for 
the usefulness of th e mu ltipulse control of spin dephasing 
are given in section IV 



II. 



A GENERAL FORMULA OF THE 
MULTIPULSE CONTROL 



We consider the pure dephasing phenomena of a \ spin 
caused by an interaction with a reservoir. Here and hence 
force, we use the term "spin" as an arbitrary two-level 
quantum system with energy difference tko. The Hamil- 
tonian of this system is 



"Hr = Ho + 7~Lsb i 

where H.q and Hsb are defined by 

Ho = hujS z + Hb 
Hsb = hS z B . 



(1) 



(2) 
(3) 



Here S z is the z-component of the ^ spin, Hb is a Hamil- 
tonian of the reservoir and B is an interaction amplitude 
operator which consists of reservoir operators and causes 
the energy change in the energy difference between two 
states due to reservoir. 

It is our purpose to study effectiveness of a multipulse 
control on this spin system. Assuming that pulse dura- 
tion is so short that the spin-reservoir interaction is ig- 
nored during pulse excitation, we write the corresponding 
Hamiltonian Hp during a pulse as 



Hp{t)=Ho + W P {t) 



(4) 



where Wp(t) denotes: 

W P (t) = -^H ■ (S+e- iuJt + S„e luJt ) . (5) 

Here we consider the case where an applied square pulse 
of the magnetic field H is on resonance with the spin 
whose magnetic moment is defined by p. Using an elec- 
tric field instead of the magnetic field H and a dielectric 



moment instead of the magnetic moment /2, we can ap- 
ply the following formulation to an electronic state in a 
quantum dot. 

The phase information of the spin is obtained by calcu- 
lating the off-diagonal element of the reduced density op- 
erator, (l|Tr^ p(t)\0) where |1)(|0)) is the upper (lower) 
state of the spin, p(t) is the density operator of the to- 
tal system and Ttr denotes the trace operation over the 
degrees of the freedom of the reservoir. When we apply 
regular multiple pulses with an interval r s and a pulse du- 
ration At after applying a pulse at an initial time (t — 0), 
the quantity is written as 

(l\Tr R P (t)\0) 

= Tr R {l\e- lLR ^- {NTs+At ^T + [e~ lS ^s + ^ dt ' Lp - ]{t ' ) ] 

N-l 

x{~[[ e^ LR ^-^T + [e- l ^^ + ^ dt ' Lp ^ t ' ) ]}p{0)\{)) , 



3=0 



(6) 



where T + is the time ordering symbol from right to left, 
and L v (y — {P, j} or {R}) is the Liouville operator 
defined as 



-\H 



v i J J 



and Hpj is the Hamiltonian during the j-th pulse. 

For later convenience, we introduce the hyperspace 
notation as |m)(n| — > \mn)), Tr[(\k)(l\y \m)(n\] — > 
{(kl\mn)) and Tr[(\k) (l|)t(0 \m)(n\)\ -> ((kl\0\mn)) 
for an arbitrary hyperoperator O. Using the complete- 
ness relation, 



E \ mn ))(( mn \ = 1 



(8) 



m=0,l n=QS 

Eq.(^|) is separated into the following elements as 

((W|e _i Jn dt ' Lp{t,) \mn)) 
= e- iw « fc - i ) i2 -( m -")* 1 >(cos(^)4, m +isin(^)(l - 4,™)) 



x (cos(-)5„,j - i sin(-)s(l - S n j)) 
with the so-called pulse area 9 = I* 1 ft • H, and 

—iLftt 



(9) 



((kl\e 1 a pR\mn)) = Sk, m 5i, n exp[-iL min t]pR . (10) 



In Eq.(|10[), pr is the density operator for the reservoir, 
and 



with 



L m . n PR = HB,m PR — PR Hb,u , (H) 



2m — 1 

Hs,m = = (hu + B)+H B - (12) 



Now we consider the case where we first apply a J 
pulse to generate a superposed spin state at an initial 
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time (t = 0) and then apply N times of tt pulses with 
constant interval t s . Assuming the initial total system 
to be 



p(p) = PH (P)\0)(0\ 
and using Eqs. ^ ~ we have for even N, 
(l\Tr R P (t)\0) 



(13) 



= Tr R [((10\ 



-iL R (t-Nr s ) 



|io» 



x{((01| e -' i ^|01))«10|e-^^|10))}^ / 9 fl 



= Tr R [t 



iL li0 (t-NT B ) | e -!Xo,i r ag — iii.o r a 



} 2 H 



5 i«B,o T Sc ,iH B ,i T s jf p in B ,o (t-N-r, 



(14) 



In the last line of Eq.(14), the exponentials with ex- 



ponent Hb,i and TCb,o alternately appear, which reflects 
the switching between states |1) and |0) by tt pulses. This 
is schematically drawn in Fig.pl The diagram consists of 
two parallel straight lines: the upper(lower) line is asso- 
ciated with forward (backward) time evolution of p R in 
Eq. (|l4|) . The arrow indicates the direction of time evo- 
lution. The thick (thin) line segment corresponds to the 
time evolution in the state |1) (|0)). 




t = 



FIG. 1: The diagram for the time evolution of the reduced 
density operator in Eq.(|l4|). 

We can rewrite Eq.(|l4|) in the form 
(l\Tr R p(t)\Q) 

= (~l) N ^ e-<" ^- N ^Tr R [U(r a , (t-(N- l)r s )) 

xC/t( , t s )UH2t s , 3t s ) • • • tft((jv - 2)t s , (N - l)r s ) 
XPr(0) 

xU((N - 3)r„ (N - 2)r.) • • • U(3t„ 4t s )[/(t s , 2t s ) 
xU^(0,r s )}, (15) 

where we use the relation 

e -l(H B +B)(U-t 2 )) = e -*Wat lt ^( tlj ^gjtWofa ( 16 ) 

and 



U(t u t 2 ) = T + [exp[-~jf dt'B(t')]}, (17) 



B{t) = e^ HBt Be 



(18) 



Using the diagram method in Appendix A to take up 
to the second order of B(t) in cumulant expansion, we 
obtain the signal intensity for the coherently oscillating 
spin described by the off-diagonal element of the reduced 
density operator, (l\Tr R p(t)\0), as 

I(t) = |(l|Tr fl ,^)|0)| 2 

N 

= exph2((-l) JV £(-l)"(4n - 2)S((N -n+ l)r.) 

n=l 

N 

+ (-l) A, - 1 ^2(-l)"+ 1 5(i-nr s ) 



where we introduce the second cumulant as 



(19) 



S(t) = f dh I' dt 2 Re[(B(ti)B(t 2 ))] . (20) 
Jo Jo 



Although the formula (|19|) is obtained up to the second 
order of B{t) in cumulant expansion, we should recall 
that the formula can explain the pure dephasing in many 
systems where spin- reservoir interaction is composed of 
a strong linear interaction and a relatively weak nonlin- 
ear one. This is because (1) the higher-order cumulants 
vanish irrespective of the coupling strength for the non- 
degenerate two-level system linearly interacting with the 
boson reservoir p9|, and (2) they are negligible for the 
nonlinearly interacting system if the interaction is suffi- 
ciently weak. 

The quantity S(t) can be rewritten as 



s(t) = I 2^H- 



cos tut 



(21) 



where J{ui) is the power spectrum defined by 

J(w) = / dte iiAjt (B(t)B(Q)) . (22) 



The obtained general formula (|19|) describes the spin 
dynamics under multipulse control in terms of the cor- 
relation function of the interaction between system and 
reservoir without specifying the form of interaction am- 
plitude B, and can be applied to various systems. 

At the last part of this section, we consider the asymp- 
totic behavior of the quantity S(t) in the long time region. 
Rewriting Eq. (pl|) into the form as 



S(t) = 



00 du> 1 — cosuit 

J ( W ) ~2 



— oo 

oo 



2tt 

du) 1 — cos tot 
~2 ntuj 2 



J(w) t 



we have an asymptotic form as 



(23) 



(24) 
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since 



1 — cos tot 

■Kti0 2 



(25) 

means 
by the 



in the long time region (t 3> r c ). The relation 
that the quantity S(t) for large t is determinec 
value J(0), the power spectrum at to = 0. Using the 
definition of -^tt^ = t^-, we find that the general formula 
( |19| ) under the multipulse control with the long interval 
is reduced to 



I(t) = \Tr R (l\p(t)\0)f 



exp[-— t] . (26) 

1 2 



The multipulse control does not work in this situation. 
However, as the pulse interval becomes short, the effec- 
tiveness of pulse control increases. We will show it by 
evaluating concrete cases in the next section. 



III. MODEL CALCULATION 

Now we apply the general formula obtained in the 
previous section to the following two cases: (1) a phe- 
nomenological two-level model where the correlation 
function of the interaction decays exponentially, and (2) 
a spin-boson model where a spin nonlinearly interacts 
with a boson reservoir. 



A. Phenomenological Two-Level Model 

In order to consider the general behavior of spin dy- 
namics and to find model-independent results, we first 
apply the expression (^) to a spin system interacting 
with the reservoir where correlation function exponen- 
tially decays as 



(B(ii)B(t 2 )) = A 2 exp[-|*i - t 2 \/T c 



(27) 



is not effective. However, the multipulse control with 
sufficiently short pulse interval works well as shown be- 
low, because the memory effect of the reservoir plays an 
important role in the short time region. Equation (|2§| ) 
gives the dependence of ln[/(t)] at the pulse applied time 
t/v(= Nr a ) for even N in the form 



ln[/(f 



JVJ 



T, 



■{t N 



1 + e" 



2N) t, 



-2e" 



2e~ 



+ 4JV t, 



1 + e T c 



1 



(30) 



In the case of large r s (3> r c ), we obtain the following 
approximate expression 



ln[/(tjv)] ~ -— t N , 

introducing an effective dephasing time: 

T 

rpe — 2 

2 "(l+(w-2) T t) 



(31) 



(32) 



Comparing Eq.(pl|) and Eq.(|26|), we find that the effec- 
tive dephasing time XJf increases with decreasing to 
give suppression by the multipulse control. 

Using the dimensionless time scale normalized by the 
dephasing time T 2 as f c(s) = t c{s) /T 2 and t N = r c ( s )/T 2 , 
we evaluate the signal intensity ln[/(£jv)] as a function 
of a pulse application time for several values of the pulse 
interval f s , while the correlation time f c is fixed to 0.02 
(see Fig.|). 



where r c denotes the correlation time of the reservoir and 
A is the strength of the interaction. In this case, Eq. (|2p| ) 
is written as 



S(t) =-A 2 T c {t-T c (l 



-t/T 



)} 



(28) 



Although the correlation function ( p7|) is simple, we 
can find various features depending on the ratio t/r c . In 
the long time region (t 3> t c ), the time evolution of the 
spin is characterized by 
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S(t) 



To 



(29) 



where we defined the dephasing time T 2 = (A 2 t c ) 1 . 
(Using Eq.(||) and Eq.@, we have J(0) = 2A 2 r c = 
■Sr.) Since the memory of the detailed information on 
the interaction with the reservoir is lost, the time evo- 
lution becomes fully irreversible and multipulse control 



FIG. 2: The signal intensity ln[/(fjv)] as a function of a pulse 
application time for several values of the pulse interval fs ■ The 
pulse interval f s is set to 0.005, 0.1, 0.5 for a constant correla- 
tion time {f c = 0.02). The dashed line shows the exponential 
decay with T2 as a reference. 
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When the pulse interval f s is much larger than the cor- 
relation time (f s = 0.5), the decay is nearly determined 
by the pure dephasing time T 2 . As f s decreases to 0.1, 
we can see that the decay of the ln[/(i/v)] becomes slow. 
This means that the pulse train recovers the phase co- 
herence. When the pulse interval becomes much shorter 
than the correlation time (f s = 0.005), we can see that 
the pure dephasing is highly suppressed. 

Evaluating numerically effective dephasing time T| 
from the decay time at sufficiently large time t = 10, 
we show the pulse-interval dependence of the effective 
dephasing time T| in FigJ|. We find T| becomes large 
for sufficiently small pulse interval and monotonously de- 
creases to an asymptotic value 0.5. 
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FIG. 4: Time evolution of the signal intensity ln[/(t)] un- 
der the multipulse control. The pulse interval f s is set to 
0.05, 0.1, 0.2 for a constant correlation time (f c = 0.02). 
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FIG. 3: The pulse- interval dependence of the effective dephas- 
ing time T 2 e . 

Let us consider in more detail the physical origin why 
the dephasing is suppressed by the pulse-train control. 
The general formula ([H]) gives the time evolution of the 
system between the pulses as well as the values on the 
pulse-applied time. Figure || shows the time dependence 
of ln[/(f)] during the multipulse control for f s = 0.05, 0.1 
and 0.2 with a constant correlation time (f c = 0.02). We 
can see the decay recovers between the pulses. This is 
because the applied ir pulses cause the time reversal for 
spin dynamics. The degree of recovery reflects to what 
extent the memory of the interaction with the reservoir 
remains, which depends on the ratio — . 

The memory effects of the reservoir in optically ex- 
cited systems have been extensively studied during the 
past couple of decades as the non-Markovian optical 
problem (27], |(| || §|, || |§ " 



3§ ^ ^ ^ |3Sg. We 

should stress that the dephasing control by a 7r-pulse 
train discussed here is one of the most typical and clear 
phenomena arising from the finitcncss of the correlation 
time. 



written in an exponential form, we have quantitatively 
evaluated ln[/(t)] and discussed the physical origin of 
the suppression of the pure dephasing. Next we con- 
sider a microscopic case where the pure dephasing oc- 
curs by the interaction between a ^ spin and a reservoir 
which consists of bosons. While the standard spin-boson 
modeljl^] has been often used to describe pure dephasing 
by linear spin-reservoir interaction, there exist many de- 
phasing phenomena which cannot be explained without 
non-linear spin-reservoir interaction. A typical example 
of this situation is the pure dephasing phenomena of a 
localized electron in a solid. The dephasing phenomena, 
which are originated from the interaction with phonons, 
cannot be described by assuming the interaction to be lin- 
ear where absorption and emission of one phonon occurs. 
As shown in sec. II, the long time behavior is described by 
the value of the spectrum J (to) at u> = (see Eq.(p4|)). 
In the case of linear interaction, the value J(0) corre- 
sponds to the interaction with the acoustic phonon at 
the wave number k = 0, which corresponds to the trans- 
lation of the whole solid without distortion. Since such 
a motion does not affect the electron, the pure dephas- 
ing cannot be described with the standard spin-boson 
modelJIH- We need to consider a nonlinear interaction 
with the reservoir. The reason why the linear spin-boson 
model cannot generally describe the pure dephasing is 
given in Appendix B. 

In this paper, assuming the reservoir to be composed 
of the infinite number of bosons 



(33) 



in an equilibrium state with temperature T, we take into 
account the quadratic interaction between the spin and 
reservoir as 



B. Nonlinear Spin-Boson Model 

Applying the formula (|l^) to the case where the corre- 
lation of the interaction between system and reservoir is 



H SB = hS z B 



= yS z ^ ^2 h lnW^l^rn(bl + b\)(b ri 



I m 



bl 



(34) 
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The interaction Hamiltonian (|34|) describes simultaneous 
absorption and/or emission of two phonons with frequen- 
cies toi and <jj m where interaction strength is hi m . The 
term b,b m in the above interaction Hamiltonian causes 
the emission of a phonon with frequency loi and the ab- 
sorption of a phonon with frequency u m at the same time. 
The energy change of this process is h(u>i — u) m ), and 
therefore it contributes to J(co — 0) as shown below. 
This is the reason why the nonlinear interaction is essen- 
tial when we consider the pure dephasing phenomena. 

The quadratic interaction gives the correlation func- 
tion in Eq.(p2|) of the form 

(S(t)S(O)) 

= E E h % lm uiu m {{n{ui) + l)(n(uj m ) + 1) e ~^'+^ t 

I 771 

+2n( W| )(T l K) + l)e^'^ i 

+n(w,)n(w m ) e^'+^j , (35) 



where n(uo) is the boson distribution function: 
n(uo) 



huj/k B T ^ 



(36) 



We rewrite the summation over the label of the boson 
in Eq.(^5|) into an integral form to give the power spec- 
trum J(uf) in Eq.(22) by using the following relation, 



l m 

= E E h ^ S ( e ~ w (5 ( e ' - ^m)/(e, e') 

i n 

< poo 

de / de'h(e)h(e')f(e,e') (37) 
Jo 

where h(io) is a spectral density of the coupling func- 
tion for the spin-reservoir interaction and f(io,Lo') is an 
arbitrary function of co and a/. 

When the spectral density of the coupling function is 
assumed to be a Gaussian distribution with mean fre- 
quency uip and variance 7 p ; 



He) 



exp(- 



(e - ujp) 2 



(38) 



we have 

JM 

= 2tt J deh(e) {e{uj - e)(n(e) + l)(n(w - e) + l)h(u) - e) 

-2e{uj - e)n(-(ui - e))(n(e) + l)h(-(u - e)) 

-e(w + e)n(e)n(-(w + e))ft(-(w + e))} . (39) 

Figure | shows J(u>) for % = 0.4, Sq = s 2 = 1/40, 
and hiOplkftT — 1 where we have defined as Co = lo/lo p 
and 7 p = ^f p /uip. We can see three peaks: The right (left) 
peak at u> = 2(= —2) in Fig.^ corresponds to the emission 



(absorption) of two bosons by the spin. The center peak 
at Co — represents the simultaneous boson absorption 
and emission, which causes the pure dephasing. This is 
because the quantity S(t) is described as S(t) ~ ^p-t in 

the long tine region, and the quantity gives the de- 
phasing rate for the dynamics (see Eq.(M) and Eq.(|26])). 




FIG. 5: Power spectrum J (Co) for j p — 0.4, Sq = 1/40, and 
hio p /k B T = 1. 

Now we discuss the effectiveness of the multipulse con- 
trol in suppression of the pure dephasing. In Fig.0, we 
show In [7(f)] for the same conditions as in the Fig.g, us- 
ing a definition as t = to p t. The decay is well suppressed 
for small pulse interval f s (= lo p t s ) = 0.1, which is consis- 
tent to the results in the previous subsection. However, 
a surprising result is obtained in increasing the pulse in- 
terval: As the interval increases to f s = 1.5, the degree of 
the suppression becomes worse, but a significant recovery 
is seen for f s = 2 and 2.5. 
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FIG. 6: Logarithmic display of intensity In [7(f)] for several 
values of pulse interval f s = 0.1, 1.5, 2, 2.5 in the case of 7 P = 
0.4, Sq = 1/40, and Tuo p /k B T = 1. 

Next we evaluate the pulse-interval dependence of the 
effective dephasing time T% which is obtained by the de- 
cay time at i = 45. Figure ^ shows that the character- 
istic feature inherent to the dynamical behavior of bo- 
son system. When the pulse interval is near to — , the 
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more effective suppression of pure dephasing is caused 
by resonance enhancement between the oscillatory be- 
havior of boson system and the pulse train. This should 
be contrasted with the phenomenological case where T| 
monotonously decreases with increasing the pulse inter- 
val t s (see Fig.||). 
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FIG. 7: The pulse interval dependence of the effective de- 
phasing time T-2 for the non-linearly interacting spin-boson 
model. 



IV. CONCLUDING REMARKS 

We have obtained a general formula (|l^) which de- 
scribes the time evolution of spin system under pulse- 
train control. The formula is expressed with the cor- 
relation function for arbitrary interaction amplitude B 
which causes dephasing. This enables us to clarify how 
the suppression of pure dephasing depends on the pulse 
interval and the characteristic time of the interaction. We 
have applied the general formula to two examples: (1) a 
phenomenological two-level model where the correlation 
function of the interaction decays exponentially, and (2) a 
spin-boson model where a spin nonlinearly interacts with 
a boson reservoir. In the case (1), the analytic form for 
the intensity is evaluated from the off-diagonal element of 
the reduced density operator at pulse applied times. We 
define effective pure dephasing time T| from this form 
and we find that T| effectively increases by decreasing 
the ratio of the pulse interval to the correlation time. 
Moreover, the analytic expression of the time evolution 
of the intensity between pulses enables us to generally 
discuss the mechanism of the suppression of dephasing 
by multipulse control. In the case (2), we also found 
the suppression of the pure dephasing. However, special 
attention should be paid to the fact that the extent of 
the suppression does not monotonously decrease as the 
pulse interval increases, which is contrasted with the phe- 
nomenological case (1). When the pulse interval is near 
to the half of the average period of boson system, pure 
dephasing is effectively suppressed by resonance between 
the oscillatory behavior of boson system and the pulse 



train. This coherence recovery shows us that the micro- 
scopic mechanism can be used for effective suppression. 
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APPENDIX A: DIAGRAM METHOD 

In order to obtain Eq. ( |l9| ) , we expand the time-ordered 
exponential in Eq.(|l5|) and take up to the second order 
of Bit). In the case of two-pulse excitation , for the most 
simple case where ir pulse is applied once at t — r s after a 
^ pulse at t = 0, the intensity I(t) is obtained as(27|, |2£| 

J(t) = cxp[-2(2S(r s ) + 2S{t - r s ) - S(t))\ . (Al) 



Using the fact that the exponent in Eq.(Al) is composed 
of the function S with three kinds of time difference, 
r s ,t — t s , and t, we can write a diagram as in FigM 



(a) o 



t = 0,' Ts 



(b) o .^r\ (°) o 
=» • — * ,» • — 

t t = o Ob,-' t t = o_ . 



1 ' 



FIG. 8: Diagrams for Eq.flAj). 
2S(t s ), (b) corresponds to 2S(t 
-S(t). 



Diagram (a) corresponds to 
— t s ), and (c) corresponds to 



In this diagram, the dots are placed at the pulse ap- 
plied times and the observation time t. The broken lines 
are connected between the time points that are the ends 
of the broad lines. The coefficient of the function S can 
be determined by the number of the broken line Nb and 
the number of the time points N t included between the 
broken line as (—l) Nt Ni,. 

For application of it pulses in two times, we obtain the 
intensity as 

I{t) = exp[-2(6S'(r s ) - 2S(t - r a ) + S(t) 

-2S(2t s ) + 2S(t - 2r s ))] . (A2) 

The diagram is written in Fig.^. 

These diagram rules enable us to obtain a general for- 
mula Eq.(|l9[) for an arbitrary number of pulses. It is 
also noted that the diagram method is a powerful tool to 
evaluate more general cases with irregular pulse interval 
and pulse area. 
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(a)-,- (b) (c) r - v 

t 1 / V 1 1 '° 1 ' ■ 1 N . 

t = V ^v v " 2^ t t = o x Sf- 2^-" t t = o Xs 2xs t 

o » i / o ' o i o " o i o " 




FIG. 9: Diagrams for Eq.(A2). Diagram (a) corresponds to 
6S{t s ), (b) -2S(t-r 3 ), and~(c) S(t), (d) -2S(2r s ) , and (e) 
2S(t-2r s ). 



Using Eq.(pT|), we obtain the second cumulant S(t) in the 
form 



S(t) = / dul(u)(2n{uj) + 1) 
Jo 



1 — cos mt 



(B5) 



This corresponds to T (t) of Eq.(17) in@, which de- 
scribes the decay of the off-diagonal element of the den- 
sity operator by the interaction with the reservoir. In the 
long time region, the quantity S(t) goes to 



APPENDIX B: NOTE ON THE CONVENTIONAL 
LINEAR SPIN-BOSON MODEL 



7r[/(0)(2n(0) + l)] 



t 

T 2 



(B6) 



The pure dephasing has been explained with the spin- 
boson model |l9| where a | spin linearly interacts with a 
reservoir which consists of bosons as, 



Use = hS z B = hS z J2(9*kh + 9kb\) 



(Bl) 



which shows us that the value J(0) deteremines the long 
time behavior. 

Now we consider the spectral density I(u>) of the reser- 
voir to be the one defined by Caldeira and Leggett [[l9| as 
in previous studies pi O, 



In Eq.(Bl), g k is the coupling strength between the spin 
and the fc-th boson of the reservoir. For this model, the 
correlation function in Eq.(j2^) is obtained as 

(B(t)B(O)) = J2 \9k\ 2 {(n(w k ) + l)e-~ fc * + n^e^ 1 } , 
k 

(B2) 

which gives the power spectrum of the form, 

J(to) = 2ir{I(Lu){n(uj) + l)d{uj) + I{-u))n{-u)9{-uj)} . 

(B3) 

Here we have introduced the spectral density of the cou- 
pling function I(uS) as 



/>oo 

V] \g k \ 2 5(uj - uj k ) ■ ■ ■ = / duil(u>)- 
i Jo 



(B4) 



(B7) 



where a determines the coupling strength with the reser- 
voir and an integer n depends on details of the system- 
reservoir interaction[|l9). Evaluating T2 with Eqs.(fB6|) 
and (B7), we have 



To 



2h 

irk B T 



(n=l) 
00 (n 7^ 1) 



(B8) 



which shows that the pure dephasing is absent except for 
the ohmic dissipation case (n = 1). It is remarked that 
the model can describe the pure dephasing only for the 
ohmic dissipation case. 
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